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The  Derivation  of  a  Multivariate 
Probability  Density  Function  Having  an 
Exponential-Type  Bivariate  Marginal  Density 


1.  BACKGROUND 

This  report  discusses  the  derivation  of  a  multivariate  probability  density 
having  certain  desired  characteristics  related  to  the  statistical  characterization 
of  terrain  features,  particularly  for  use  in  calculation  of  the  scattering  of  elec¬ 
tromagnetic  waves  from  uneven  terrain  surfaces.  A  brief  discussion  of  some 
surface  scattering  aspects  will  serve  to  place  this  report  in  its  over-all  context. 
The  initial  effort  in  this  area  is  described  in  some  detail  in  a  report  by  Lennon 
and  Papa.  1  The  application  of  the  technique  to  the  actual  data  base  outlined  in 

that  report  and  the  subsequent  use  of  the  output  in  electromagnetic  scattering 

2 

calculations  have  resulted  in  reasonable  agreement  with  data.  The  probabilistic 

3 

elements  enter  in  the  representations  of  terrain  cross  sections.  Ruck  et  al 
give  expressions  for  the  average  bistatic  rough  surface  cross  section  oq  which 


(Received  for  publication  2  May  1980) 

1.  Lennon,  J.F.,  and  Papa,  R.J.  (1980)  Statistical  Characterization  of  Rough 

Terrain,  RADC-TR-80-9.  - 

2.  Papa,  R.J.,  and  Lennon,  J.F.  (1980)  Statistical  analysis  of  digitized  terrain 

maps  and  its  application  to  predicting  cross  section  of  rough  terrain, 
Proceedings  of  the  Second  Workshop  on  Terrain  and  Sea  Scatter.  George 
Washington  University,  Washington,  D.  C. 

3.  Ruck,  G.T.,  Barrick,  D.C.,  Stuart,  W.D.,  and  Krichbaum,  C.K.  (1970) 

Radar  Cross  Section  Handbook,  Vol.  2,  Plenum  Press,  New  York. 
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depend  on  the  assumed  two-point  probability  density  function  (PDF)  relating 
heights  of  involved  surface  elements: 


The  height  distributions  affect  J,  the  factor  associated  with  the  contribution  from 
specularly  oriented  surface  slopes.  For  the  case  where  the  surface  heights 
(Zj,  z,j)  have  a  bivariate  Gaussian  relationship. 


Alternatively,  if  the  heights  are  assumed  to  share  a  bivariate  exponential  form 
PDF, 


P,^‘  (;*>)  •'3l5Vs,,/2 


then 


J  = 


/  2  -2\1/2 

/  3T2  \ 

('/6T\  (?x  +  Cy\ 

(  2 .2  )  exp 

V  2ct  )  \  .2  f 

V®  5z  > 

\  Cz  / 

where 


a  -  variance  of  the  bivariate  PDF, 

z  =  (  1  \  =  two  dimensional  height  vector, 

\  z2  / 


covariance  matrix  for  the  heights, 
correlation  length  for  the  surface, 
sin  -  sin  9 
0.0, 

cos  -  1.0  . 
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The  second  term  of  this  cross  section  relation  consists  of  the  scattering  matrix 
elements  3^.  This  becomes 

(1  +  cos  2a)  R..(a) 

*w  s  IcoFTT^osVl  Vertical  polarization) 
i  a 


(1  +  cos  2a)  Rjte) 
(cos  +  cos  flg) 


(horizontal  polarization)  . 


In  these  expressions. 


cos  <r  -  v  e  -  sin  a 

rmc>  =  - 7^=f=; 

«r  cos  a  +  V  er  -  sin  a 


cos  «  -  Ver  -  sin  a 


cos  a  +  Vtf  -  sin  a 


where 


0 ^  =  angle  of  incidence  (with  respect  to  surface  normal), 
=  angle  of  scatter  (with  respect  to  surface  normal). 


•■(W  • 


Here,  «r  is  the  relative  complex  dielectric  constant  of  the  surface;  the  subscript 
II  refers  to  the  E-fleld  in  the  plane  of  incidence;  and  the  subscript  j.  refers  to  the 
E-fleld  normal  to  the  plane  of  incidence.  These  simplified  forms  of  Ruck's 
expressions  follow  from  the  assumption  that  the  receiver  is  far  from  the  trans- 

A 

mitter  so  that  the  portion  of  the  "glistening  surface"  which  contributes  to  the 
diffuse  multipath  is  a  long,  narrow  strip  extending  between  the  transmitter  and 


4.  Beckmann,  P.,  and  Spizzichino,  A.  (1963)  The  Scattering  of  Electromagnetic 
Waves  from  Rough  Surfaces.  Macmillan  Co. ,  Mew  Yorlc. 


receiver.  This  assumption  allows  us  to  make  the  approximation  that  the  azi¬ 
muthal  scattering  angle  4g  =  0.0,  which  leads  to  the  simplified  forms. 

The  preceding  discussion  related  to  a  Gaussian  two -point -height  PDF  which 
is  the  ordinary  case,  and  an  exponential  form  found  to  give  good  agreement  for 
some  cases.  Both  situations,  however,  require  the  use  of  statistical  forms 
describing  the  relation  between  heights  on  a  two-point  basis.  In  the  initial  attempt 
to  describe  terrain  height  relations,  *  the  emphasis  was  on  consideration  of  a 
multivariate  relationship  for  the  surface  heights.  Each  region  of  interest  was 
divided  into  a  grid  structure  with  ten  points  in  each  dimension;  the  observed 
heights  at  these  grid  points  were  then  considered  as  a  single  multivariate  obser¬ 
vation  of  an  one  hundred  variate  joint  PDF.  The  variates  are  assumed  to  have 
equal  means  and  variances.  Determination  of  the  type  of  multivariate  PDF  which 
is  more  appropriate  for  the  observed  data  set  was  made  by  a  hypothesis  test. 

For  the  Gaussian  case,  there  is  no  problem.  It  is  well-known  that  the  form 
of  the  Gaussian  is  preserved,  whatever  the  order  of  the  multivariate  PDF.  Thus 
the  bivariate  marginal  density  for  an  N -variate  Gaussian  (here  N  =  100)  is  always 
a  bivariate  Gaussian.  However,  it  is  not  clear  what  form  would  be  required  as 
the  multivariate  PDF  when  an  exponential  form  is  the  bivariate  marginal  density. 
Unlike  the  Gaussian,  the  form  of  the  PDF  for  that  case  is  not  preserved  for  suc¬ 
cessive  marginal  densities.  This  requires  some  additional  explanation. 

In  that  basic  report,  1  it  is  shown  that  an  N-variate  exponential  PDF  would 
have  the  form 


P<7', 


ZN>  = 


(N+l) 


N/2 


N+l 
,  2 


<2») 


N-l 

2 


1/2 


-[(N+1KzTR"1z)]1/2 


Note  that  this  contains  the  gamma  function,  T(x)  ■  (x  -  1)1  It  can  also  be  shown 
that  the  form  for  the  general  L-variate  marginal  density  for  this  N-variate  PDF 
is  (see  Appendix  A  of  this  report  for  the  derivation) 
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Here,  K^iu)  is  the  modified  Bessel  function  of  the  second  kind,  of  order  u.  Thus 
an  N-variate  exponential  PDF  corresponds  to  a  bivariate  marginal  density 


2 


and  not  to  a  simple  bivariate  exponential  PDF.  Since  the  form  of  the  N-variate 
PDF  that  would  have  the  exponential  bivariate  marginal  density  was  not  immed¬ 
iately  evident,  the  initial  assumption  was  that  the  surface  heights  were  either 
multivariate  Gaussian  or  exponential  in  distribution.  The  question  of  the  form 
for  a  multivariate  PDF  which  would  reduce  to  the  bivariate  exponential  marginal 
density  was  left  for  future  consideration. 

In  this  report,  the  derivation  of  the  general  N-variate  PDF  having  a  bivariate 
exponential  marginal  density  is  presented. 


2.  APPROACH 


This  section  will  describe  the  basic  approach  used  to  establish  a  general 
form  for  a  multivariate  PDF  having  a  bivariate  exponential  marginal  density. 

The  results  are  based  on  extrapolation  from  lower  to  successively  higher  order 
variates  until  a  trend  is  evident.  In  the  following  section,  the  hypothesized  form 
for  the  N-variate  PDF  will  be  verified. 

Prior  to  the  specific  analysis,  some  comments  of  a  general  nature  should  be 
made.  In  the  original  report,  *  it  was  shown  that,  for  convenience,  the  correlated 
form  of  the  quadratic  form  appearing  in  the  PDF  can  be  replaced  by  an  uncorre¬ 
lated  form  through  a  linear  transformation  of  the  variates.  A  further  simplifica¬ 
tion  in  variates  was  introduced  which  reduces  the  form  to  one  having  unit  variance. 
This  procedure  will  be  adopted  here,  leading  to  the  equivalent  relations 


P^z •  •  •  >  = 

p(yr  •  •  •  <  =  ci 


”^c2^zl’ zn^  5  • 

^"fc2^1yl+*2y2+*  ’”+*'NyI^ 


■fc2^wl+w2+*  *  *  +WN^ 


P(w, 


V  =  C1  • 


1/2 


’  *  *  ZN^ 

1/2 


1/2 
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The  results  are  the  same  for  the  nonzero  mean  case.  The  V  values  are  the 
eigenvalues  of  the  inverse  of  the  covariance  matrix.  The  value  of  the  quadratic 
form  is  unchanged  by  these  transformations.  In  the  determination  of  the  appro¬ 
priate  functional  forms  (as  distinct  from  the  actual  calculations  of  PDF  results) 

2 

the  uncorrelated  (o  =  1)  form  will  be  considered  directly. 

One  additional  point  which  should  be  discussed  is  that,  in  the  above  relations 
and  in  the  remainder  of  this  report,  the  covariance  matrix  R  and  consequently 
its  inverse  are  positive  definite  symmetric.  This  is  a  requirement  of  the  prop¬ 
erties  of  probability  densities.  For  the  present  case,  it  is  satisfied  by  our 
assumption  that  the  correlation  _£  has  a  Gaussian  spatial  dependence  and  that  the 

elements  of  the  corresponding  covariance  matrix,  a..  c(R..l  have  the  form 
2  i  ‘3 

=  p  .^.<7  .  In  the  report  by  Lennon  and  Papa  it  was  shown  that  when  the  ele¬ 
ments  have  that  form,  the  conditions  for  R  to  be  positive  definite  are  satisfied. 
Once  R  is  positive  definite,  then  so  is  R"*.  This  result  assures  convergence 
for  the  distribution  integrals. 

The  initial  extrapolation  from  a  bivariate  exponential  form  was  the  result 
of  manipulations  associated  with  determining  coefficients  for  the  N-variate 
exponential  PDF.  *  This  involved  transformation  of  coordinates  to  an  (r,  6)  space. 
Consider  the  function: 


f(Wj,  w2,  w3,  w^)  =  Cj(Wj  +  w2  +  w^  +  w2)-1/2  e  2 


\/w2+w2+w2+w2' 


This  is  reduced  to  a  function  of  two  variables  by  integration: 


f(wj,w2)  =  Cj  j  j  <»;  +  ...  +w2> 


-1/2  -Cg's/ w2+. . .  +w^' 


dwg  dw^ 


Now  let 


2  2  2  .  2  2  2 
r  =  w3  +  w^  and  a  =  w  ^  +  w2 


Then 


«  -c2(r2+a2)1/2 

f(wl,w2)  =  cx  J  f  - 


J  J  /  5  o  N 

0  0  •»  +  a^ 


dr  d« 
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and  setting 


2  2  * */2 

du  =  (r  +  a  )  r  dr 


gives 


This  shows  that  a  function  of  the  assumed  type  reduces  after  integration  to  an 
exponential  function  in  two  variates.  Since  the  major  interest  in  the  first  terms 
is  the  progression  to  more  complex  forms,  the  determination  of  the  coefficients 
required  to  produce  a  PDF  from  the  general  functional  form  at  each  level  will  be 
discussed  in  Appendix  B. 

The  form  that  reduces  to  the  bivariate  exponential  is  a  combination  of  a 
rational  function  and  an  exponential  term.  This  four-variate  relation  suggests 
the  assumption  that  subsequent  forms  also  consist  of  combinations  of  such  terms. 

The  six-variate  form  is  assumed  to  be 


Then 


or  with  the  same  change  in  variables 


■  •!  /  de/[(^7)F 


■C'  vW' 


where 


2  2  ,  2 
r  =  w5  +  w6 


l2  =  W^  +  Wg  +  Wg  +  W^  . 


If  we  now  let  u  =  Vr2  +  a2  we  have 


f'<wr  w2,  w3,  w4)  =  2*^  J  CjU 


o  “C„u  f  o  _C2U 

~2  e  2  du  +  J  cju  1  du 


f'(w1(w2,  w3,w4)  =  2Jrc 


rc'lc2  [-Ei  (-c2'/=?)]  ; 


f'(w1(w2,w3,w4)  = 


/  2  ,  2  ,  2  2 
'W1  +  w2  +  w3  +  w4 


-c2^/w2+w2+w2+w2' 


This  verifies  that  the  assumed  form  is  appropriate.  These  integral  relations  and 
those  used  below  are  all  found  in  Gradshteyn  and  Ryzhlk.. 

Intermediate  order  forms  (P  =  3,  5)  can  now  be  determined  by  simple  single 
integrations  of  the  corresponding  higher  ones  (  p  =  4,  6).  As  before,  consider: 


«W3>-  /  (  7=  t=t=?) 

-*  \ Vwl  +  w2  +  w3  +  w4  / 


/2  2  2 

"c2  v  w  l+w2+w3+w4 


dw4  . 


5.  Gradshteyn,  I.S.,  and  Ryzhik,  I.M.  (1965)  Tables  of  Integrals,  Series,  and 
Products,  Academic  Press,  New  York. 


2  2  2  2  2  2 
Then,  with  r  =  w^  +  a2  and  a  =  wi  +  w2  +  w3 


f(Wj,  w2,  Wg)  =  2 


-e2r 

•  *  dr 


'l*W2*W3)  =  2cl  Ko  (C2V^ 


2  2  2' 
1  +  W2  +  W3 


where  K  is  the  modified  Bessel  function  (MBF)  of  order  zero, 
o 

Before  addressing  the  case  of  p  =  5,  we  first  note  that  f(w w2>  Wg,  w^)  is 
also  readily  expressible  as  an  MBF: 


f(wj,  w2,  w3,  w4)  = 


_ 2cl c2 _ 

»  +  w2  +  W3  +  W4 


X  Ki/2(c2V^i  +  w2Vw2  +  w5') 


and  after  some  manipulation,  so  is  f(Wj,  ....  wg): 

c 

f(wi*  w2*  w3*  w4’  w5’  we)  =  7 — 7==nwr 

L(c2VwJ  +  ...  +  w2  j 

X  K3/2  (c2‘n/*i  +  •••  +  '"q) 

(Note,  cQ  is  the  constant  resulting  from  the  changed  form. )  Then 

/*  T  co 

H»1.w2.w3.w4.w5)  =  J  ,  :  — gv\3/2 

-»  L(c2V«'1  +  ...  +w6  j 

X  Kg/2  (c^w2  +  ...  +  w2'  )  dw6 


2  2  2  2  2 
The  substitutions  a  =  +  ...  +  and  x  =  (Wg/a  +  1)  give 


f(w 


l,w2.  w3,  w4,w5>  =  [(a2)'1/2  <c^f3/2  co]  f  (x‘1/2)  /  (x-l)‘l/2 

X  K3/2  (c2  ^ 


and 


flwj,  w2,  w3,  w4,  Wg)  = 


— k.  (ci 

-2  /TT  1  \ 2V  1  5 

'2  vwl  +  •••  +  w5  / 


At  this  point,  functional  forms  have  been  constructed  starting  with  the 
bivariate  exponential  and  progressing  through  stages  to  that  for  the  six -variate 
case.  If  specific  forms  for  the  corresponding  PDF  relations  (as  evaluated  in 
Appendix  B)  are  now  introduced,  an  extrapolation  to  the  general  term  of  the  pro¬ 
gression  can  be  estimated.  Note  that  the  correlated  form  of  the  variates  is  now 
used: 


p*zl>  z2) 

/  3  \ 

\2»le2l1/2/ 

• 

pfZj,  Zg,  Zg)  * 

(  33/2  \ 

\2»2  IS3ll/2/ 

Ko  (^VD1/a) 

p(zl>  z2’  z3-  z4*  = 

(  32 
V23'2  .5'2 la. 

jT72^  (^uTsilfJl/2) 

XK,/;!(^li 

W/2) 

dx 
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ptej,  z2,  z3,  z4,  zg>  = 


-3/2 


X  Kt  (^/TlzTg*l^l 1/2  ) 

/  3^  \  /  t  »i  \ 

p(Zj,  z2#  z3,  z4,  z5,  z6)  =  ^  25/2  #7/2  jl/2  )  ^  ^6  / 

XK3/2(^[zTfi;lzi1/2)  . 

Based  on  these  forms  the  hypothesis  is  advanced  that  the  general  form  is: 
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It  now  remains  to  be  shown  that  this  satisfies  the  criteria  for  an  N-variate  PDF, 
and  that  it  has  the  desired  bivariate  exponential  marginal  density. 


3.  MULTIVARIATE  PROBABILITY  DENSITY  FUNCTION 

In  this  section,  the  general  multivariate  form  will  be  shown  to  satisfy  require¬ 
ments  for  a  probability  density  function  (PDF).  As  discussed,  the  covariance 

T  -1 

matrix  is  positive  definite  symmetric  and  hence  the  quadratic  form  [z  g  zj  is 
always  positive;  so  is  the  function.  It  remains  to  be  shown  that  the  integral  of 
this  function  is  unity. 

This  involves  evaluating  the  integral: 

oo  00 

lo  *  f  •  •  •  J  P(zj>  •  •  • »  zj^)  ^zi  •  •  •  *^zn 


i 


This  result  confirms  that  the  general  multivariate  form  actually  is  a  valid  repre¬ 
sentation  of  a  PDF. 

The  demonstration  that  this  PDF  does  have  a  bivariate  marginal  density  can 
be  shown  directly.  However,  an  approach  which  leads  to  a  more  general  result 
will  be  employed  here.  The  expression  for  an  arbitrary-order  marginal  density 
will  be  derived  and  subsequently  evaluated  for  the  bivariate  marginal  case. 

The  L-variate  marginal  is  given  by 


p(z1(...,z  )  = 


/•••  /  /-ErsT - WTrUil/2>'(¥) 

-co  -oo  ^2  2  »  2  |r|1/2/ 


«„  (^uW/2)  dz^  . . .  dz 


L+l  * 


After  coordinate  transformation  and  the  substitution 
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i  -  wN  +  w^_ j  +  . . .  +  w^+j  and  a  =  w  ^  +  Wg  + 


,.  +  w 


p(w.,  . . . ,  w.  )  = 


\  /  l  ""-“I 
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X  Kn_3  ^/3  a2')  ds  . 
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The  introduction  of  a  normalized  variable  X  =  (s  /a  +1)  leads  to 
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N-L  1 

2 

L 

l  N-l  N+l  1 

\  J,T  / 

V 

N-2L-3 

dL 

N-3 


ce  N-L-l  . 

/ '*  - 11  (i)  k¥ 


19 


After  simplification,  this  becomes 
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Then 
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and  finally. 
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dx  . 


This  is  the  expression  for  the  L-variate  marginal  density  of  the  given  PDF. 
To  confirm  that  it  has  the  desired  bivariate  form,  let  L  =  2.  Then, 


^UT£2'V1/2 


UTR 


When  the  relations  K_^(x)  =  K^(x)  and  =  (*/ 2x)*/^  *'x  are  introduced, 

this  becomes 


To-l.i1/2 


/ 3  7*2  \  “^UTE21*) 

‘i-H'-y—w-3)  • 


If  we  take  |r2I^2  =  Xg)  then  we  can  see  that  this  is  the  desired  bivariate 
marginal  form,  and  the  general  multivariate  PDF  does  indeed  satisfy  all  the 
requirements. 

4  CONCLUSION 

The  multivariate  probability  density  function  (PDF)  which  would  have  a 
bivariate  exponential  marginal  density  has  now  been  determined.  At  this  point, 
the  digitized  terrain  data  can  be  analyzed  with  the  new  multivariate  form,  replac¬ 
ing  the  original  multivariate  exponential  PDF  in  the  hypothesis  testing  procedure. 

In  addition,  this  PDF  can  be  employed  in  alternative  future  terrain  characterization 
techniques.  Finally,  there  still  remains  the  question  as  to  just  how  much  the  elec¬ 
tromagnetic  relations  for  terrain  scattering  are  affected  by  the  use  of  this  form, 
as  opposed  to  the  multivariate  exponential  PDF. 
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Appendix  A 


Determination  of  the  L-Veriate  Marginal  Density  of  a 
Multivariate  Exponential  Density 


The  purpose  of  this  section  is  to  determine  the  general  expression  which 
specifies  the  marginal  density  of  any  order  for  an  N-variate  exponential  probabil¬ 
ity  density  function.  The  procedures  are  similar  to  those  of  the  Appendices  A 
and  B  of  the  report  by  Lennon  and  Papa.  * 

The  multivariate  correlated  exponential  PDF  is: 


p(z. ,  ...,  = 


(N  +  1)N^2 _ 

Ebl 


uW/2 


To  obtain  the  L-varlate  marginal  we  integrate  over  (N  -  L)  variates: 
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This  represents  coordinate  transformations  to  an  uncorrelated  set  of  variates 
2 

with  or  =  1.  The  integration  next  uses  the  Jacobian  relation  to  transform  to  an 

2  2°  2 

(r,  ....  space.  Let  r  =  wN  +  wn-1  +  +  wl+1;  then 
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and  we  have 
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For  the  purposes  of  integration  we  have 
2*  x 
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Appendix  B 

The  Transformation  of  Functional  Forms  into 
Probability  Densities 


The  purpose  of  this  section  is  to  determine  values  for  the  coefficients  that 
will  satisfy  the  requirements  for  a  probability  density  function  for  the  functional 
forms  with  successively  two  through  six  variates.  The  procedures  are  similar 
to  those  used  in  the  report  by  Lennon  and  Papa1  involving  the  zeroth  and  second 
moment  integrals.  In  addition,  the  comments  on  R  and  the  modified  Bessel 
function  properties  also  apply  in  general  and  will  not  be  discussed  for  each 
particular  case.  For  each  case,  R^  and  its  associated  eigenvalues  (X)  are  the 
appropriate  ones  for  that  order. 

The  bivariate  exponential  case  is  known  as 


p(zlf  z2)  = 


(2*  ls,l1/z) 


The  case  for  N  =  3  starts  from: 


Rzj,  *2*  ^3^ 
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Then  as  before,  the  first  moment  is  equal  to  unity  or 
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-00  -  00  -oo 


or,  after  transformation. 
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Similarly,  for  the  second  moment 
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These  integrals  lead  to, 
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respectively,  or 
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For  N  =  4,  the  form  is 
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The  two  moment  integrals  here  are 
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The  corresponding  moment  integrals  are 


The  evaluation  of  these  integrals  leads  to 


The  moment  integrals  are 
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After  the  same  type  of  evaluation  sequence,  the  results  are 


b  =  >/3  and 
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This  completes  the  sequence  of  successively  higher  order  variate  terms  used  to 
extrapolate  to  the  general  N-variate  form  of  the  PDF. 
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